Metastable anisotropy orientation of nematic quantum Hall fluids 
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We analyze the experimental observation of metastable anisotropy resistance orientation at half 
filled quantum Hall fluids by means of a model of a quantum nematic liquid in an explicit symmetry 
breaking potential. We interpret the observed "rotation" of the anisotropy axis as a process of 
nucleation of nematic domains and compute the nucleation rate within this model. By comparing 
with experiment, we are able to predict the critical radius of nematic bubbles, R c ~ 2.6/j,m. Each 
domain contains about 10 4 electrons. 
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I. INTRODUCTION 

Quantum liquid crystals are gapless condensates that 
spontaneously break rotational and/or translational sym- 
metry. There is by now a large amount of theoreti- 
cal and experimental work studying these new phases of 
strongly correlated systems in different realizations like 
quantum hall systems^, high T c superconductors^ and 
heavy fermion compounds^. 

The two dimensional quantum smectio 4 -, also referred 
to as stripe phase, is a metallic state that breaks transla- 
tion invariance in one direction. It was conjectured that 
this modulated electronic configuration is a good ground 
state of two dimensional electron gases (2DEG) un- 
der specific values of an external magnetic field^&i&ii. 
In fact, at partial filling factor the system tends to 
separate into homogeneous fluids with different densi- 
ties. Coulomb repulsion frustrates this tendency and 
the system is forced to rearrange itself by lowering its 
dimensionalit y 10 ' 11 . Collective excitations of this ground 
state were computed and anisotropic gapless fermionic 
correlations, very different from those in the usual the- 
ory of Fermi liquid s 10 : 11 : 12 : 13 : 14 , were obtained. 

Strong thermal, as well as quantum fluctuations of 
the stripes, could produce topological defects, disloca- 
tions or disclinations that, under appropriate circum- 
stances, can melt the stripe order into an homogeneous 
but anisotropic liquid 15 . This new state is called the 
quantum nematic stat o 16 : 17 , and it is probably the best 
candidat e) 18 : 19 to explain the anisotropics observed in 
2DEG at half filled Landau levels&£. In fact, experiments 
are compatible with the interpretation of a spontaneous 
rotational symmetry breaking at approximately l50mK 
and a weak native potential, responsible for aligning 
the principal axis resistance, of the order of ImK per 
electron^ -. On the other hand, no pinning was detected 
in the I — V curves, which suggests a liquid state rather 
than periodical arrays. 

Some models 2 ^ were proposed to understand the native 
symmetry breaking potential responsible for the align- 
ment of the anisotropy. However, the origin of this po- 



tential remains unknown. In a recent experiment^, the 
structure of the native potential was studied on highly 
mobility samples over a large scale of temperatures and 
magnetic fields. It was reported a non trivial behav- 
ior of the resistance anisotropy as a function of tem- 
perature and filling factor. It was also found that the 
"easy" direction (the direction with lower resistivity) can 
be aligned along the (110) or the (110) crystallographic 
axes of the host GaAs structure. The actual direction 
preferred by the system depends on the filling factor and 
on the in-plane magnetic field. These directions can be 
interchanged according to the magnetic field sweep in 
such a way that an interesting hysteresis pattern, typical 
of metastability, comes up close to half filling. Moreover, 
to confirm the picture of a bi-stable potential, Cooper 
et. alS^ were able to quench the system in a metastable 
state (say the "easy" direction along (110)). Then, they 
observed the slow relaxation to the equilibrium state, 
aligned with the axis (110), for several final tempera- 
tures. 

In this work, we analyze this result in the framework 
of a nematic quantum fluid, submitted to a two com- 
ponent external potential; one component with nematic 
symmetry and the other one with tetragonal symmetry, 
possibly induced by the host GaAs structure 2 -. For this 
purpose, we introduce an XY model, describing an ef- 
fective quantum nematic phase, with a general external 
symmetry breaking potential that produces a two orthog- 
onal minima structure. Within this picture, we assume 
that the decay of the metastable state can be produced 
by thermal activation over a potential barrier. We ex- 
pect that bubbles nucleation of the true ground state 
into the metastable state, produced by long wavelength 
thermal fluctuations, is the main mechanism responsible 
for the decay. This assumption is reasonable provided 
the energy of the critical bubble is much greater than 
the equilibrium temperature. Wc will show that this is, 
in fact, the case for the data of ref. [U 

We use the Langer— homogeneous nucleation theory to 
compute decay rates of the metastable state. We calcu- 
late the critical bubble profile of the model, and estimate 
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the critical energy and radius by using two methods: an 
analytical variational approach and an exact numerical 
computation. By comparing with the experiment, we are 
able to predict the radius of critical domains of the order 
of 2.6pm, containing approximately 10 4 electrons. We 
find that the time evolution of the anisotropy resistance 
observed in ref. [22| is in agreement with the picture of 
nucleation of nematic domains, with directors pointing 
along a stable direction, in a metastable nematic back- 
ground with the principal axis aligned in the perpendicu- 
lar direction. We also show that the homogeneous nucle- 
ation theory is good enough to make estimations at first 
order, while the thin wall approximation is not quite ac- 
curate. Our calculations predict that the domains have 
broad walls in the actual experimental conditions. 

We present our model of quantum nematic in an exter- 
nal potential in section [TTJ To make the paper self con- 
tained and to fix notations, we briefly review the theory 
of two-dimensional homogeneous nucleation in our con- 
text in 31111 Then, in SjTV] we compute the decay rate, 
the critical energy and the radius of the critical bubble 
as a function of the parameters of our model, by using 
a variational approach and the thin wall approximation. 
To check our approximations, we numerically integrate 
the differential equation which defines the critical bub- 
ble and compare the results with our previous analytical 
estimations in section [V] Finally, we compare our anal- 
ysis with experimental results in WII and summarize our 
conclusions in section I VIII 



II. XY MODEL OF A QUANTUM NEMATIC 
LIQUID IN A SYMMETRY BREAKING 
POTENTIAL 

A nematic state is a homogeneous orientational or- 
dered state with the forward and backward directions 
identified. That means that if the system has a prefer- 
ence axis orientated along an angle 9, the state has the 
nematic symmetry 9 — > 9 + it. In two dimensions, this 
property is encoded in the definition of the complex order 
parameter Q = p e l26 , where the argument 29 guarantees 
the nematic symmetry. Close to the isotropic-nematic 
transition, we can expand the free energy in powers of Q, 
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where a 2 and 04 are arbitrary constants which depend 
on the microscopic details of the system. V(hQ) is an 
explicit symmetry breaking term which depends on some 
external field h. Differently from three dimensions, in 
the two dimensional Ginzburg-Landau expansion of the 
nematic order parameter there is no cubic term. 

In the absence of external symmetry breaking poten- 
tial, the modulus of the order parameter spontaneously 



gets a non zero value p for a 2 < 0, while the angle 9 
remains arbitrary. In two dimensions, the angle fluctu- 
ations are logarithmically divergent^. Therefore, in the 
absence of external fields there is no true order, but alge- 
braically decay quasi-long range order. However, in the 
presence of a small symmetry breaking potential, this di- 
vergence is removed. In this case, even small values of 
the external potential can produce big values of the order 
parameter due to the huge susceptibility of this transi- 
tion. 

The dynamics of the lowest energy modes is governed 
by the coarse-grained Hamiltonian (as usual, we have as- 
sumed at low temperatures a constant modulus of the 
order parameter in eq. (|2.ip ) 
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where J is the typical energy scale of the Kosterlitz- 
Thouless^ phase transition, and V(9) is an arbitrary 
potential that explicitly breaks rotation invariance, but 
preserves the nematic symmetry in such a way that 
V(9) =V(6 + n). 

It is possible to parametrize this potential in terms of 
its even Fourier coefficients 



V(6) = h 2n cos(2n9) 



(2.3) 



where h 2 is related with a nematic external field, /14 is a 
symmetry breaking coefficient with tetragonal symmetry 
and so on. 

This model, with n = 1 in cq. (|2.3[) . was used to 
fit the isotropic/ anisotropic transition of a Hall liquid 
at v = 9/2 filling factor, by means of Monte Carlo 
simulations 1 ^. It was shown that the general picture of 
a nematic liquid in a small symmetry breaking field cor- 
rectly describes the transition. To use this model for 
fitting experimental data, it is necessary to relate the ne- 
matic order parameter with observables. In ref. [lH it 
was shown that the relation 
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is a good approximation over a huge temperature range, 
except at extremely low temperatures where quantum 
fluctuations become important. In eq. (|2.4p . p xx and p yy 
are the measured longitudinal resistivities in the x and y 
direction respectively. 

In this paper, we adopt the same criteria and study 
higher harmonics of the external field. For this purpose, 
we will analyze the effect of the second term n — 2 in 
cq. (12. 3|) . Thus, we will consider a potential of the form 



V{9) = h 2 cos(26») - h 4 cos(46») , 



(2.5) 



where h 2 > and /14 > are coefficients that measure 
the relative weights of the nematic and tetragonal com- 
ponents of the host symmetry breaking potential. It is 
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FIG. 1: Potential of eq. (|2.3|) with hi < 4/i 4 . The metastable 
states are at 6 = nir, while the true ground state is align 
with 9 = 7r(2n + l)/2 with n = 0, ±1 . . .. Maxima are taken 
at cos(20Max) = h^/iAhA,). The energy difference between 
minima is 2/i2 and the height of the barrier in related with hi 
through eq. (|2.6p . 



interesting to note that, while the term hi is somehow 
mysterious, the term h 4 is not prohibited by symmetry 
and could be induced by the square structure of the host 
GaAtM. 

The net effect of the new term h 4 is to introduce local 
metastable minima for the angular variable 9. We depict 
the potential for typical values of hi < Ah^ in fig. [TJ 
The sign of /14 controls which one of the different minima 
is metastable and which one is the true ground state. 
We arbitrarily choose this sign in such a way of having 
a metastable state at 9 — mr and the ground state at 
9 = n(2n + l)/2 with n = 0,±1 . . .. On the other hand, 
the potential has maxima at cos(20Max) = hi/(4h4,). 

The energy difference between the stable and 
metastable minima is 2h 2 , while /14 is related with the 
height of the energy barrier, in fact, 
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which in the limit of quasi degenerate minima reduces to 
V(8 Max ) - V(0) = 2/14. 

With this potential it would be possible, by a quench- 
ing process, to prepare the system in a metastable phase; 
for instance, a nematic state with principal axis pointing 
in the 9 — direction. In these conditions, we would ex- 
pect that long wavelength thermal fluctuations could pro- 
duce nucleation of nematic domains with principal axis 
in the perpendicular direction 9 = Tt/2. Therefore, we 
should see that the anisotropy of the resistivity "rotates" 
or changes in time from one direction to its perpendicular 
one as it has been observed^. 



III. TWO DIMENSIONAL BUBBLE 
NUCLEATION 

One of the possible mechanisms for thermal activated 
decaying is the nucleation^ of ground state domains in 
a homogeneous metastable state which fills all the avail- 
able area. The dynamics is determined by the domain's 
energy that, in general, is a competition between a bulk 
contribution (proportional to its area) and a boundary 
term (proportional to its perimeter). 

Consider, for instance, a bubble with spherical sym- 
metry of radius R as depicted in fig. [5J In the thin wall 
approximation, that is, when the width of the wall is 
much smaller than the radius, the boundary and bulk 
contributions to the energy are well defined, 
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where AT is the energy difference between the stable and 
metastable states per unit area, a is the surface tension, 
and the ellipsis indicates subleading order in the thin wall 
approximation. 

While for small radius the positive boundary term 
dominates the energy, for large R, the negative bulk con- 
tribution dominates. Due to this competition, there is a 
critical radius R c — a j AT, where the energy has a max- 
imum E c — -no 1 j AT, and the critical bubble is at un- 
stable equilibrium. The supercritical bubbles (R > R c ) 
will grow until filling all the area with the ground state. 
On the other hand, the subcritical bubbles (R < R c ) will 
shrink and finally disappear. Both contributions are im- 
portant in a phase transition since the actual mechanism 
is given by random long wavelength thermal fluctuations 
which generate all types of bubbles. Some of them will 
grow and others will shrink. In this picture, the transi- 
tion is completed when the true ground state percolates 
the metastable one. The relative importance of these 
contributions depends on the probability of fluctuations 
and on the growth rate of the supercritical bubbles 
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The important quantity to study nucleation is the nu- 
cleation rate per unit of area. In the homogeneous nu- 
cleation theory of Langer 24 , this quantity is given by 
r = VCD exp(—E c /T), where E c is the energy of a critical 
bubble, T is the final equilibrium temperature, f2 is the 
growth rate of a slightly supercritical bubble 2 ^ and the 
prefactor T> comes from the computation of fluctuations 
around the critical bubble profile. 

The computation of Q and T> from microscopic quan- 
tum models is a very difficult task (see ref. H3 and refer- 
ences therein). However, it was shown 2 ^ that in two di- 
mensions and in the thin wall approximation (AT — * 0) 
these quantities can be cast in terms of the macroscopic 
parameter AT, 
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FIG. 2: Sketch of a symmetrical bubble which contains a 
nematic liquid with the director pointing in the 9 = n/ 2 di- 
rection. It is embedded in a metastable state composed of a 
nematic liquid with the director pointing to 9 = 0. p is the 
modulus of the nematic order parameter and n is the director. 



Interestingly, due to its two-dimensional character, this 
expression has no corrections in powers of the thin wall 
adimensional parameter— ATT /a 2 . 

Homogeneous nucleation theory is reliable provided 
E c » T. For E c ~ T, small amplitude thermal fluc- 
tuations could trigger the phase transition and nucle- 
ation and spinodal decomposition can no longer be dis- 
tinguished. 

Provided the conditions for homogeneous nucleation 
are satisfied, the probability of nucleating several bubbles 
at the same time is very small. We can estimate the 
typical time to complete a transition as the time of a 
simple nucleation event. Then, 
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where r is the time estimated to complete the transition, 
r is the nucleation rate per unit of area, and A is the 
total area of the sample considered. 

This "static" approach will be sufficient for the purpose 
of this paper. However, time dependent corrections to 
Langer theory can be evaluated using out of equilibrium 
Shwinger-Keldish techniques^. 

In the next section we will compute the time r as well 
as the coefficients AT and a in terms of the parameters 
of our model, J, ft 2 and /14. 



IV. NEMATIC CRITICAL BUBBLES AND THE 
THIN WALL APPROXIMATION 

A critical bubble is a radially symmetric static field 
configuration that solves the following differential equa- 
tion 

1 dV{6) 
J 80 

with the boundary conditions linv—,00 9{r) = and 
lim^oo 9'(r) — 0. V(9) is the potential of eq. 
(ng.P). 
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FIG. 3: Typical radial symmetric bubble profile that solves 
eq. 14.11 The continuous line represents the bubble configura- 
tion which starts close to the true ground state 9 = tt/2 and 
reaches the metastable state 9 = at asymptotically large ra- 
dial distance. The dash line is the square derivative (properly 
rescaled) that defines the critical radius and the wall thickness 



Such a solution describes a "bubble" like configura- 
tion which starts close to the true ground state 9 = tt/2 
and reaches the metastable state 9 = at asymptoti- 
cally large distances. The change from the stable to the 
metastable state occurs around the critical radius R c , 
over a distance £ which defines the wall thickness of the 
bubble. We depict a typical profile of a critical bubble 
and we also draw the square derivative of the profile in 
fig. [31 The maximum of the derivative defines the criti- 
cal radius and the width of the peak is a measure of the 
bubble wall thickness that, as we will show, is related to 
the nematic correlation length in the metastable phase. 

Eq. (|4.ip is an extremely difficult differential equation 
to solve analytically and we will show a numerical treat- 
ment in the next section. However, it is possible to have 
some insight of its behavior through a variational anal- 
ysis. The idea is to propose a reasonable ansatz for the 
solution, by considering the critical radius and the wall 
thickness as variational parameters. Then, we determine 
these parameters by extremizing the critical energy. We 
make the following ansatz, 



1 — tanh 



Rr 



(4.2) 



where the radius R c and the wall thickness £ will be de- 
termined extremizing the energy. 

The form of this function is inspired in the problem of 
an asymmetric quartic potential. In that case, eq. (|4.2p 
is the exact solution of eq. (|4.ip for the one-dimensional 
problem and, for higher dimensions, it is the correct form 
in the thin wall approximation £/R c << 1—. Although 
our potential is much more complicated than a simple 
quartic one, we expect to grasp the general behavior 
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FIG. 4: Comparison between a numerical solution of eq. (I4.1|) 
(continuous line) and the variational ansatz eq. (|4.2|l (dash 
line) . 



and the correct order of magnitude with this ansatz. Of 
course, a variational technique is not a well controlled ap- 
proximation, thus, in the next section, we compare our 
variational analytical results with a numerical treatment 
of eq. (|4.ip . Just to have some feeling on the approxima- 
tion, we compare the variational profile of eq. (I4.2| with 
the exact one, obtained from direct numerical integration 
of eq. gjjj, in fig.H 

Substituting eq. (|4.2[) into the Hamiltonian eq. (|2.2p . 
we have the energy of the bubble as a function of the 
variational parameters, 

E(R c ,0= I d 2 x J - V6 b (r,R c ,$ 2 + V(9 b ((r,R c ,£)))- 



(4.3) 

Integrating to leading order in £/ R c we find an expression 
similar to eq. (|3.1[) . 

E(R c ,0 = -nATR 2 c + 27Ta(C)R c + 0((aRc) 2 ) , (4.4) 
where 



AT = V(0) - V(tt/2) = 2h 2 



(4.5) 



essentially comes from the integral of the potential. On 
the other hand, ct(£) has contributions from both terms 
of the integral and is given by 
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The first term comes from the gradient contribution and 
grows when the thickness narrows. On the contrary, the 
second term, which comes from the potential, is a linear 
increasing function of £. Therefore, we fix this parameter 
looking for a stationary solution 
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and obtain the optimal value, 
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which is two times the correlation length on the 
metastable phase. With this value for the thickness, the 
superficial tension gets the simplest form, 



~ = J^ 2 ^/Ah 4 -h 2 . 



(4.9) 



Now, we can determine the critical radius R c by im- 
posing 



dE(R c ,Z) 
dR c 

We have, in this way, 



2ttR c AF + 2ttct(C) = 0. (4.10) 
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Finally, plugging AjF, a and R c into eq. (|4.4[) . we obtain 
the energy of the critical bubble, 



E c 7T J 4/l 4 

T = 2 \~h7 ~ 1 



(4.12) 



It is necessary to have in mind that the thin wall ap- 
proximation also impose restrictions on the values of h 2 
and h^, since 

This completes our estimation of relevant dynamical 
quantities in terms of the parameters of our model in the 
thin wall approximation. 



V. NUMERICAL CALCULATIONS 

In order to check the range of validity of our variational 
approach we numerically evaluate the bubble profile and 
the energy. 

Assuming a solution with radial symmetry and rescal- 
ing the variables with the magnetic length l c as 
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L 2 A — » U 2 ,4 



2 2 n2 

r — > r irt c , 



(5.1) 



we find the adimensional equation 
d 2 9 1 dO 



dr 2 r dr 



2h 2 sin(2f5) - Ah A sin(46») = , (5.2) 



<•) = and 



with the boundary conditions linv^r 
linWx, 0'(r) = 0. 

To solve this equation we transform the boundary 
value problem into an initial condition problem and use 
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the shooting method to seek for solutions. In this method 
one shoots the initial derivative, integrates the equation 
with a Runge-Kutta method of order 4 and, according 
to the accuracy to match the boundary values, the ini- 
tial conditions are corrected. This procedure is iterated 
looking for convergence. 

To actually solve the equation, we need to fix hi and 
hi. There is by now extended experimental work in dif- 
ferent samples and regions of magnetic field and den- 
sity that clearly shows that, while the anisotropy appears 
around T = 150mK, the energy scale of the aligning po- 
tential is about lmK per electron. In our model, this 
is compatible with values of hi and hi of the order of 
1CP 2 in units of J/n^- For instance, in ref. [IH, the 
isotropic/anisotropic transition at v = 9/2 was success- 
fully fitted with a value of hi — 0.05. Of course, the 
specific value may changed for different samples and for 
different filling factors and in-plane magnetic fields. As 
we have stated before, due to the bidimcnsionality, the 
prefactor in the decay rate, eq. (|3.3p . depends just on 
the energy difference between the minima, say 2/i2, and 
not on the height of the barrier, proportional to /14. This 
means that for a given value of the decay time, the energy 
of the critical bubble only has a logarithmic dependence 
on hi. Then, what actually matters for the calculation 
of the energy of the critical bubble is the order of hi 
and not its precise value. For this reason, we will fix 
a reasonable value of hi and we will make the numeri- 
cal calculations for several values of /14 in the range of 
10~ 2 . When comparing with experiments, we will com- 
ment again on these values and will show the robustness 
of the results for small changes of these quantities. 

In fig. [5]we show a set of solutions with fixed hi = 0.02 
for several values of /14 while in fig. [6] the square modu- 
lus of the derivatives are depicted for the same values of 
the parameters. Firstly, we observe the general features 
of a critical bubble described in the previous sections as 
expected. All solutions begin near 9 = it/ 2 and asymp- 
totically go to 9 = 0. The critical radius R c can be read 
from the position of the maxima of fig. [6] and we can es- 
timate the wall thickness as the width of peaks at some 
arbitrary fixed height (half the peak height, for instance). 
According to the previous section, we see that R c is an 
increasing function of /14 and, the greater /14 the smaller 
the wall thickness. We collect this information in fig. [Jj 
where we draw the critical radius as a function of /14. The 
continuous line connecting the points is a polynomial fit. 
We also draw with a dash line the function R c (h4) given 
by eq. (|4.1ip . Interestingly, except for very low values of 
/14, the variational result reasonable agrees with the exact 
result; the greater /14 the better the approximation. 

Finally, we compute the energy of each critical bubble 
by numerically integrating eq. (|4.3| . In fig. [8j we show 
the energy as a function of the critical radius. The con- 
tinuous line is a polynomial fit of second order. Note 
that we are perfectly fitting five points with a second 
order polynomial. The reason for that is simple: the 
critical energy can be cast in terms of the critical radius 




FIG. 5: Bubble profiles for different values of hi. We have 
fixed h2 = 0.02 for all solutions, r is measured in units of 
•Jvlc and /14 in units of J/-kH. 2 c . 



as E c = AT 7ri? 2 . We have estimated this parameter 
to be AT — 2h 2 . Therefore, we expect a quadratic de- 
pendence of the form E c — a 2 R^,, with the quadratic 
coefficient given by a 2 = 2hiir. We have found a fitting 
value ai — 0.1257, in excellent agreement with the thin 
wall estimation with the value hi = 0.02, fixed for all the 
critical bubbles. It seems striking to compare figures [7J 
and [5] While in fig. [5] the exact results match thin wall 
calculations almost perfectly, in fig. [7J we see a clear de- 
viation. The reason is that in determining AT ', the only 
contribution comes from the second term of eq. (|4.3p . 
The contribution of order R\ essentially comes from the 
constant part of the bubble profile and our variational 
function estimates quite well this area. Indeed, this is the 
only source of error in the results in fig. [5] and, for this 
reason, the numerical calculations and the variational es- 
timations match almost perfectly. On the other hand, 
for the estimation of a, necessary to evaluate the critical 
radius (see eq. (|4.11|) ). we have two contributions (see 
eq. (|4.6[) ). Both of them come from the bubble wall and, 
as we can see from fig. [4j our variational ansatz is not 
so good in that region. However, as it is shown in fig. [JJ 
the approximation gets better as the value of ft.4 grows, 
since the barrier height increases and, consequently, the 
wall gets thincr in this limit. 



VI. EXPERIMENT INTERPRETATION 

In this section we analyze the experimental results of 
ref. I22I in the context of our model of nematic nucleation. 
The main purpose of this section is to analyze whether 
the above mention data can be interpreted as thermal 
activation over a barrier and if the homogeneous nucle- 
ation theory of Langer is applicable. Moreover, we want 
to check the idea of a nematic liquid proposed earlier—^ 
to describe the ground state properties of these systems. 
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FIG. 6: Wall profiles for different values of /14. We have fixed 
fi2 = 0.02 for all solutions, r is measured in units of \pnl c 
and /14 in units of J 




FIG. 7: Critical radius as a function of /14. The continuous 
line is a polynomial fit. The dash line is the critical radius in 
the thin wall approximation calculated with eq. (I4.11[l , We 
have fixed /12 = 0.02. R c is measured in units of ^/n£ c and fi4 
in units of J/-k£ 2 c . 

To reach this aim, we need to make not only qualitative 
comparisons with experiment, but concrete predictions, 
specially about the size of the nematic domains and walls, 
in order to further check this picture. 

In ref. H2, Cooper et. al. showed clear evidences 
of metastable behavior in the resistance anisotropy ori- 
entation of very clean 2DEG. The measurements were 
done in a square sample of area 25mm 2 and density 
N s = 3 x 10 11 cm~ 2 . The main result was that the 
"hard" and "easy" directions of the longitudinal resis- 
tance depend, at half filling, v = 9/2, 11/2, 13/2 . . ., on 
the magnetic field sweep. A typical hysteresis diagram 
associated with metastability was shown and, in the same 
work, the authors were be able to "quench" the state into 
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FIG. 8: Critical energy in units of J as a function of the 
critical radius R c expressed in units of •Jtx(. c - The continuous 
line is a polynomial fit of order two. 

a metastable direction and to follow the slow relaxation 
to the equilibrium. The decay rate strongly depended on 
temperature and on filling factor. 

The procedure for analyzing the experiment in the con- 
text of our model is the following: we take three mea- 
sured quantities from the experiment, the time needed 
to complete the anisotropy "rotation" , the equilibrium 
temperature at which the decay was observed and the 
critical temperature for the isotropic/anisotropic transi- 
tion. With this input, we compute the energy of the 
critical bubble by inverting eq. (|3.4p . Note that this cal- 
culation is model independent and its result gives us in- 
formation about the applicability of Langer theory. With 
this value of the Energy we estimate the critical radius 
R c from fig. [5] The value of /12, as was explained before, 
was taken from several previous experiments and theo- 
retical fittings, and the value of the new parameter /14 is 
predicted form fig. [7] 

In table Q] we have collected the experimental results 
for v = 13/2 in the first two columns, r is an estimation 
of the typical time to complete the transition and T is 
the equilibrium temperature at which the time evolution 
was witness. 

Assuming that the main mechanism is thermal activa- 
tion, it is immediate to estimate, from this data, the value 
of the critical energy by using eq. (|3.4|) . The area was 
taken from the experiment and we have fixed the value 
of h% — 0.02 J I (ttI 2 .). Due to the logarithmic dependence, 
the critical energy is not sensible to this particular value 
but just to its order. We have fixed it by considering that 
the strength of the native symmetry breaking potential 
is about lmK per electron 2 -. The results are depicted in 
the third column of table U 

The first observation is that E c /T takes values between 
40 and 50. Therefore, the homogeneous nucleation theory 
of Langer is a reasonable approximation, at least for first 
order estimations since E c /T >> 1. Note that this value 
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of the energy corresponds to E c ~ 3K. 

The next step is to compute E c /J. The stiff- 
ness J can be considered of the same order of the 
isotropic/anisotropic transition, J = a T c , where a is 
a constant of order one and T c is the critical tempera- 
ture for the isotropic/anisotropic transition. While T c is 
an experimental data, a should be computed from the 
model. The order of T c ~ 150mK is not difficult to 
obtain; a more accurate number is generally more in- 
volved since the transitions are rounded by disorder and 
other effects. On the other hand, the value of a for our 
model is not simple to compute. When ignoring the small 
symmetry breaking potential, the renormalization group 
analysis of the Kosterlitz-Thouless transition^ gives an 
estimation of a ~ 3. Moreover, Monte-Carlo simulations 
of the full models, fitting the experimental data of ref. 8 
(at v — 9/2), give a reduced value of a ~ 1.1. For the 
estimation of E c / J we have used T c ~ 150mK, which is 
reasonable for v = 13/2 (note in ref. |22| that at lOOmK 
the anisotropy is completely developed at this filling fac- 
tor), and the value of a obtained from a renormalization 
group analysis. These values are shown in the fourth 
column of table [B 

Now we are ready to predict the value of the critical 
radius and the parameter /14 by using the numerical cal- 
culations of figures [5] and [7] respectively. We show these 
results in the last two columns of table HI 

We find for the critical radius R c ~ 2.6/im, where 

o 

we have considered the magnetic length £ c = 197 A (at 
v = 13/2). By taking into account the electronic density 
of the sample, this value gives an estimation of 10 4 elec- 
trons inside the critical bubble. We see that this predic- 
tion is completely reasonable since, the dimension of the 
critical domains is big enough to approximately contain 
ten broken stripes in each domain, provided we consider 
the Hartree-Fock value of the stripe period 5 as a reason- 
able estimation. Moreover, the domains are 10 6 times 
smaller than the size of the sample. 

Although the precise value of h% was fixed by hand, the 
size of the critical domain predicted is not very sensible 
to that value. In fact, since the critical energy is only 
logarithmically dependent on hi and E c = ttATR^., we 
see from eq. 14.51 that R c ~ Xj^fh-i- That means that 
when varying h% over a range of reasonable values, the 
critical radius only changes up to 15%. We also note 
that the predicted value of /14 is almost equal to h 2 and 
shows a very tiny temperature dependence. Recently, 
it was pointed out^ that piezoelectricity in GaAl can 
induce an aligning potential of the form cos(40). In that 
work, it was predicted that the barrier between the two 
degenerated minima is 10 times the Coulomb energy, 
roughly lmK per electron. This value is in complete 
agreement with the values of hi predicted in this paper. 

We have repeated the calculations with different values 
of h% and we have found a linear relation h^jhi ~ 0(1). 
This is compatible with eq. (|4.12|) which, for a fixed 
value of E c , predicts the linear dependence in the thin 
wall approximation (up to logarithmic corrections). Of 



T [mK] 


T [s] 


E c /T 


E c /J 




hi 


50 


3.6 x 10 4 


51 


5.7 


6.8 


0.020 


70 


6 x 10* 


47 


7.3 


7.6 


0.023 


90 


1 


41 


8.2 


8.1 


0.026 



TABLE I: Summary of experimental results and theoretical 
predictions. T is the equilibrium temperature of the final 
state, r is the typical time that takes the transition to be 
compieted. E c is the energy of the critical bubbie, R c is the 
critical radius, £ c the magnetic iength and hi the tetragonal 
component of the native potential measured in units of J/ir£c . 

course, in order to determine a precise value for the pa- 
rameter /12 it is necessary to fit the data for all values of 
temperatures above the critical temperature. 

It is important to note that, although the homogeneous 
nucleation approximation seems to be a reasonable one, 
the thin wall approximation is not accurate for the regime 
of this experiment. This fact can be observed in fig. [7] 
By using eq. (|4.13[) we estimate £/R c ~ 2/3. Then, the 
domain walls are not sharply defined getting broader on 
a smooth interphase of approximately 1/xm. 

VII. SUMMARY AND DISCUSSIONS 

We have analyzed the experimental observation^ 
that the native potential, responsible for the resistance 
anisotropy alignment in quantum Hall fluids at half fill- 
ing, has a non-trivial structure, which favors two orthog- 
onal directions. 

In this framework, we have studied an XY model, 
which describes a nematic fluid in an external symmetry 
breaking potential, compatible with the nematic symme- 
try. We have considered the tetragonal coefficient of the 
Fourier expansion of the potential and we have shown 
how it can produce two orthogonal local minima struc- 
ture. 

In this picture, we have assumed that the "rotation" , 
observed in the anisotropy axis, is mainly driven by ther- 
mal activation over a barrier. Then, the decay of a 
metastable direction into the orthogonal direction, which 
represents the true ground state, is dominated by nucle- 
ation of nematic domains. 

To compute decay rates we have used the homogeneous 
nucleation theory^ 4 - and, after comparing with experi- 
ment, we have concluded that it is a reasonable approx- 
imation for first order estimations since E c /T ~ 50 for 
the data of ref. I22I 

We have implemented an analytical variational ap- 
proach inspired in the bubble solutions of a quartic poten- 
tial and we have computed the critical energy and radius 
in the thin wall approximation. In order to check the 
quality and the range of applicability of our approxima- 
tions, we have numerically integrated the critical bubble 
differential equation, and we have computed its energy 
for a wide range of parameters. 

By comparing with the experiment, we were able to 
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predict the radius of the nematic critical domains of the 
order of R c ~ 2.6fim; each domain approximately con- 
tains 10 4 electrons. However, we found that the width 
of the wall is quite broad. There is a smooth transition 
from the true ground state to the metastable state spread 
in a region of the order of £ ~ 1/ito. Therefore, although 
the mechanism of homogeneous nucleation seems to be 
a reasonable assumption, the thin wall approximation is 
not accurate for the actual regime showed in the data. 

The prediction on the size of the critical domains is 
robust in the sense that it does not strongly depend on 
the detailed value of hi- On the other hand, it is quite 
sensible to the experimental determination of the critical 
temperature T c for the isotropic/anisotropic phase tran- 
sition and to the computation of the stiffness of the XY 
model in the presence of an external field. A Monte-Carlo 
fitting of the complete data, including high temperatures 
around T c for different filling factors, would help to im- 
prove the accuracy of the predictions of R c and hi- 

Finally, we would like to point out that, while the 
mechanism proposed in this paper is in agreement with 
experiment, it is not the only possibility for its explana- 
tion. Indeed, the origin of metastability could be pro- 
duced from competition between different crystal liquid 



phases rather than from a native potential or, perhaps, 
from combinations of both effects. In a recent paper—, 
it was shown that the competition between nematic and 
hexatic order parameter leads to a first order phase tran- 
sition with a clear signature of metastability in the orien- 
tation of the principal axis. The quantum hexatic phase 
could be produced from melting of the crystal bubble 
phase. All these anisotropic states are in fact competing 
close to half filling, contributing to the interesting com- 
plex structures of longitudinal resistivity reported. We 
hope to report on this possibility soon. 
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